Pseudo-Dirac Neutrinos, a Challenge for Neutrino Telescopes 
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Neutrinos may be pseudo-Dirac states, such that each generation is actually composed of two 
maximally-mixed Majorana neutrinos separated by a tiny mass difference. The usual active neutrino 
oscillation phenomenology would be unaltered if the pseudo-Dirac splittings are Sm'^ < 10"^^ eV^; m 
addition, neutrinoless double beta decay would be highly suppressed. However, it may be possible 
to distinguish pseudo-Dirac from Dirac neutrinos using high-energy astrophysical neutrinos. By 
measuring flavor ratios as a function of L/E, mass-squared differences down to ~ 10~^* eV^ 
can be reached. We comment on the possibility of probing cosmological parameters with neutrinos. 



PACS numbers: 95.85.Ry, 96.40.Tv, 14.60.Pq 

Are neutrinos Dirac or Majorana fermions? Despite 
the enormous strides made in neutrino physics over the 
last few years, this most fundamental and difhcult of 
questions remains unanswered. The observation of neu- 
trinoless double beta decay would unambiguously signal 
Majorana mass terms and hence Icpton number violation. 
If no neutrinoless double beta decay signal is seen, it may 
be tempting to conclude that neutrinos are Dirac parti- 
cles, particularly if there is independent evidence from 
tritium beta decay or cosmology for significant neutrino 
masses. However, Majorana mass terms may still exist, 
though their effects would be hidden from most exper- 
iments. Observations with neutrino telescopes may be 
the only way to reveal their existence. 

The generic mass matrix in the (i^l, {i^r)'~^) basis is 



rriL mo 



(1) 



A Dirac neutrino corresponds to the case where rriL — 
niR = 0, and may be thought of as the limit of two 
degenerate Majorana neutrinos with opposite CP parity. 
Alternatively, we may form a pseudo-Dirac neutrino 
by the addition of tiny Majorana mass terms m^, <^ 
mo , which have the effect of splitting the Dirac neutrino 
into a pair of almost degenerate Majorana neutrinos, each 
with mass ~ mo- The mixing angle between the active 
and sterile states is very close to maximal, tan(26') = 
2mo/{mR — mi,) ^ 1, and the mass-squared difference 
is (5m^ ~ 2mu{mL + mji). For three generations, the 
mass spectrum is shown in Fig.^ The mirror model can 
produce a very similar mass spectrum 0, 01 ■ 

The current theoretical prejudice is for the right- 
handed Majorana mass term to be very large, to^ ^ m^, 
giving rise to the see-saw mechanism. Then the right- 
handed states are effectively hidden from low energy phe- 
nomenology, since their mixing with the active states is 
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FIG. 1: The neutrino mass spectrum, showing the usual solar 
and atmospheric mass differences, as well as the pseudo-Dirac 
splittings in each generation (though shown as equal, we as- 
sume they are independent). The active and sterile compo- 
nents of each pseudo-Dirac pair are Vja and Vjs, and are max- 
imal mixtures of the mass eigenstates v'^ and vj . Neither the 
ordering of the active neutrino hierarchy, nor the signs of the 
pseudo-Dirac splittings, has any effect on our discussion. 



suppressed through tiny mixing angles. This is desirable, 
since no direct evidence for right-handed (sterile) states 
has been observed (we treat both solar and atmospheric 
neutrinos as active-active transitions, and do not attempt 
to explain the LSND |^ anomaly) . If right-handed neu- 
trinos exist, where else can they hide? An alternative to 
the see-saw mechanism is pseudo-Dirac neutrinos. Here, 
although the mixing between active and sterile states is 
maximal, such neutrinos will, in most cases, be indistin- 
guishable from Dirac neutrinos, as very few experiments 
can probe very tiny mass-squared differences. 

In the Standard Model, mo arises from the conven- 
tional Yukawa couplings and hence its scale is compara- 
ble to other fermion masses. In the see-saw model, 
is identified with some large GUT or intermediate scale 
mass, and thus small neutrino masses are achieved. For 
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pseudo-Dirac masses, on the other hand, we need both 
rriL and mu to be small compared to mu. The small- 
ness of with respect to m^i follows from their SU (2) ^ 
properties; the former breaks it while the latter is invari- 
ant under it. A similar property with respect to a SU{2)fi 
(obtained with a low-energy SU{2)l (S> SU{2)r symme- 
try group) may also make ma small compared to mu. 
Specific examples which achieve precisely this are given 
in Ref. . While there still remains the problem of keep- 
ing m£) itself small enough, so that the physical neutrino 
masses are tiny compared to the other fermions, there are 
a number of suggestions of how this may arise 0, 0, 0] ■ 

Astronomical-scale baselines {L > E/5m?) will be 
required to uncover the oscillation effects of very tiny 
(5m^ 0,^3- Crocker, Melia, and Volkas have considered 
possible distortions to the spectrum (ijj. Fig.[21shows 
the range of neutrino mass-squared differences that can 
be probed with different classes of experiments. Present 
limits on pseudo-Dirac splittings arise from the solar and 
atmospheric neutrino measurements. Splittings of less 
than about 10^^^ eV'^ (for vi and V2) have no effect on 
the solar neutrino flux :4| i while a pseudo-Dirac splitting 
of v-i could be as large as about 10"'* eV^ before affecting 
the atmospheric neutrinos. 

Note that models with light sterile neutrinos often con- 
flict with big bang nucleosynthesis limits on the number 
of light degrees of freedom in thermal equilibrium in the 
early universe. However, the sterile component of each 
Pseudo-Dirac pair will not be populated, provided the 
mass splitting of each pair is sufficiently small, as will be 
the case for the examples we consider here. 

Formalism. — Let (:/j^ , i/^ , t';^ ; z/j~ , i/^ , i'^ ) denote 
the six mass eigenstates, where z/+ and i'~ are a nearly- 
degenerate pair. A 6 X 6 mixing matrix rotates the mass 
basis into the flavor basis (z/g, i^^^, i^r ; J^e; ^r)- g^H" 
eral, for six Majorana neutrinos, there would be flfteen 
rotation angles and flfteen phases. However, for pseudo- 
Dirac neutrinos, Kobayashi and Lim have given an 
elegant proof that the 6x6 matrix Vkl takes the very 
simple form (to lowest order in 5m? /ra^): 





Ur 



V2 



-iV2 



(2) 



where the 3x3 matrix U is just the usual mixing matrix 
determined by the atmospheric and solar observations, 
the 3x3 matrix Ur is an unknown unitary matrix, and Vi 
and V2 are the diagonal matrices Vi = diag(l, 1, 1)/V2, 
and V2 = diag(e-*'^^e-*'^^e-*'^3)/^, The are arbi- 
trary phases. As a result, the three active neutrino states 
are described in terms of the six mass eigenstates as: 



(3) 



The nontrivial matrices Ur and V2 are not accessible to 
active flavor measurements. The flavor conversion prob- 
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FIG. 2: The ranges of distance and energy covered in various 
neutrino experiments. The diagonal lines indicate the mass- 
squared differences (in eV^) that can be probed with vacuum 
oscillations; at a given I//J5, larger &rv? values can be probed 
by averaged oscillations. The shaded regions display the sen- 
sitivity of solar, atmospheric, reactor, supernova (SN), short- 
baseline (SBL), long-baseline (LBL), LSND '3| and extensive 
air shower (EAS) experiments. We focus on the KM3 region, 
which describes the parameter space that would be accessible 
to a 1-km^ scale neutrino telescope, given sufficient flux. Cur- 
rent neutrino flux estimates for extragalactic sources indicate 
that it will be a challenge for km-scale experiments to make a 
sensitive test of the scenario proposed here, and larger scale 
experiments would likely be necessary. 



ability can thus be expressed as 



3 



(4) 

The flavor-conserving probability is also given by this for- 
mula, with j3 = a. Hence, in the description of the three 
active neutrinos, the only new parameters beyond the 
usual three angles and one phase are the three pseudo- 
Dirac mass differences, (5m| = (m^)^ — {mj)'^. In the 
limit that the Sm^ are negligible, the oscillation formu- 
las reduce to the standard ones and there is no way to 
discern the pseudo-Dirac nature of the neutrinos. 

We assume that the neutrinos oscillate in vacuum. The 
matter potential from relic neutrinos can affect the as- 
trophysical neutrino oscillation probabilities, but only if 
the neutrino-antineutrino asymmetry of the background 
is large, of order 1 For present limits on that asym- 
metry, of order 0.1 Il3l. or for less extreme redshifts than 
assumed in Ref. [l^l, matter effects are negligible. 

Supernova neutrinos from distances exceeding 
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parsecs will arrive as a 
50/50 mixture of active and sterile neutrinos due to 
vacuum oscillations. However, we focus on the poten- 
tially cleaner signature of flavor ratios of high-energy 
astrophysical neutrinos. 

_L/i?-Dependent Flavor Ratios. — Given the enor- 
mous pathlength between astrophysical neutrino sources 
and Earth, the phases due to the relatively large solar 
and atmospheric mass-squared differences will average 
out (or equivalently, decohere) . The neutrino density ma- 
trix p is then mixed with respect to the three usual mass 
states but coherent between the two components of each 
pseudo-Dirac pair: 

P = \ll^o.jZ\^.j?{\yt){4\ + (5) 

a j = l 



Here Wa is the relative flux of at the source, such that 
X^a ''^a ~ ^- "^^^ probability for a neutrino telescope to 
measure flavor Vf3 is then = (yp\p\vfi)^ which becomes 



E 



If// 



/3jl 



5r7i| L 
4E 



(6) 

In the limit that Sm'j 0, Eq. © reproduces the stan- 
dard expressions. The new oscillation terms are negligi- 
ble until E/L becomes as small as the tiny pseudo-Dirac 
mass-squared splittings Srn^. 

Since |C/e3p ~ 0, the mixing matrix U for three active 
neutrinos is well approximated by the product of two 
rotations, described by the "solar angle" 6'soiar and the 
"atmospheric angle" 0atm — 45°. The pion production 
and decay chain at the source produces expected fluxes 
of We = 1/3 and Wf_i = 2/3. In the absence of pseudo- 
Dirac splittings, it is well known 0| that this results in 
P/3 ~ 1/3 for all flavors, thus the detected flavor ratios are 
I'e : : i^r = 1 : 1 : 1. Here and elsewhere, this — i^r 
symmetry is obtained when 6'atm = 45° and Ue3 = 0. If 
pseudo-Dirac splittings are present, we thus expect 

SPp = -^[|C//31PX1 + |C/^2|'X2 + |C//33PX3] , (7) 

where SPp = — i, and we have defined, for shorthand. 



Xj 



Smj L 
AE 



(8) 



In the absence of pseudo-Dirac terms, flavor democracy 
is expected. However, the pseudo-Dirac splittings lead 
to an oscillatory, flavor-dependent, reduction in flux, al- 
lowing us to test the possible pseudo-Dirac nature of the 



TABLE I: Flavor ratios Ve ■ for various scenarios. The 
numbers j under the arrows denote the pseudo-Dirac split- 
tings, Srrij, which become accessible as L/E increases. Oscil- 
lation averaging is assumed after each transition j. We have 

used Satm = 45°, 6lsolar = 30°, aud Ue3 = 0. 



1 : 1 

1 : 1 

1 : 1 

1 : 1 

1 : 1 

1 : 1 



4/3 : 1 
2/3 : 1 

14/13 : 1 
2/3 : 1 
4/3 : 1 

14/13 : 1 



2,3 
1,3 



14/9 : 1 
2/3 : 1 
14/9 : 1 
10/11 : 1 
10/11 : 1 
2/3 : 1 



1 : 1 
1 : 1 
1 : 1 
1 : 1 
1 : 1 
1 : 1 



neutrinos with neutrino telescopes. The signatures are 
flavor ratios which depend on astronomically large L/E. 

As a representative value, we take ^soiar = 30°. Then 
the flavors deviate from the democratic -1 value by 



<5P. = -- 



SP.=SPr = -- 



(9) 



The latter equality is due to the Vp, — Vr symmetry. 

We show in Table how the : Vp, ratio is altered if 
we cross the threshold for one, two, or all three of the 
pseudo-Dirac oscillations. The flavor ratios deviate from 
1 : 1 when one or two of the pseudo-Dirac oscillation 
modes is accessible. In the ultimate limit where L/E is 
so large that all three oscillating factors have averaged 
to i, the flavor ratios return to 1 : 1, with only a net 
suppression of the measurable flux, by a factor of 1/2. 

It was recently pointed out that neutrino flavor ratios 
will deviate significantly from 1:1:1 if one or two of the 
active neutrino mass-eigenstates decay 0. The decay 
scenario bears some resemblance to that presented here. 
In particular, if there is a range of L/E values where 
the one or two heavier mass states have oscillated with 
their pseudo-Dirac partners, but the light state has not, 
then half of the heavy states will have disappeared, to 
be compared with the complete disappearance expected 
from unstable neutrinos 15j . The effects of pseudo-Dirac 
mass differences are much milder and will require more 
accurate flavor measurements than for deca vs Il5l In 
addition, the active-active mixing angles [l2| will need 
to be known independently. A detailed analysis of the 
prospects for measuring flavor ratios in km-scale neu- 
trino telescopes has been performed in Ref. jlQj. This 
study shows that it will be very challenging for km-scale 
experiments to sensitively test the pseudo-Dirac scenario, 
and larger experiments are likely to be necessary. 

Neutrinoless Double Beta Decay. — Since the two 
mass eigenstates in each pseudo-Dirac pair have opposite 
CP parity, no observable neutrinoless double beta decay 
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rate is expected. The effective mass for neutrinoless dou- 
ble beta decay experiments is given by 



2m, 



2 ^ 

3 



(10) 



which is unmeasurably smaU, (m)eff ^ 10^"' eV for the 
inverted hierarchy and even less for the normal hierarchy. 
In contrast, in the mirror model 0], the sum above has 
(m^ + mj^ , and can thus produce an observable signal. 

Cosmology with Neutrinos. — It is fascinating that 
non-averaged oscillation phases, S(j)j — Snijt/Ap, and 
hence the factors Xjj ^re rich in cosmological informa- 
tion Integrating the phase backwards in propaga- 
tion time, with the momentum blue-shifted, one obtains 



dz 



dt 



5m] H^^ 



dz 4po(l + z) 

duj 



(11) 



4po 



where Ze is the red-shift of the emitting source, and H^^ 
is the Hubble time, known to 10% This result holds 
for a flat universe, where ^Im + f^A = 1, with U,m and JIa 
the matter and vacuum energy densities in units of the 
critical density. The integral is the fraction of the Hub- 
ble time available for neutrino transit. For the presently 
preferred values flm = 0.3 and f^A — 0.7, the asymp- 
totic (ze oo) value of the integral is 0.53. This limit 
is approached rapidly: at Ze = 1 (2) the integral is 77% 
(91%) saturated. For cosmologically distant (Zf. > 1) 
sources such as gamma-ray bursts, non-averaged oscilla- 
tion data would, in principle, allow one to deduce 5m? to 
about 20%, without even knowing the source red-shifts. 
Known values of SI™ and rJA might allow one to infer the 
source redshifts Zg, or vice- versa. 

Such a scenario would be the first measurement of a 
cosmological parameter with particles other than pho- 
tons. An advantage of measuring cosmological param- 
eters with neutrinos is the fact that flavor mixing is a 
microscopic phenomena and hence presumably free of 
ambiguities such as source evolution or standard candle 
assumptions 0, 0| . Another method of measuring cos- 
mological parameters with neutrinos is given in Ref . [lOl . 

Conclusions. — Neutrino telescope measurements of 
neutrino flavor ratios may achieve a sensitivity to mass- 
squared differences as small as 10~^^ eV^. This can be 
used to probe possible tiny pseduo-Dirac splittings of 
each generation, and thus reveal Majorana mass terms 
(and lepton number violation) not discernable via any 
other means. 

Note added: As this work was being finalized, a paper 
appeared which addresses some of the issues herein l2l|. 

Acknowledgments. — We thank Kev Abazajian for 
discussions. J.F.B. and N.F.B. were supported by Fer- 
milab (under DOE contract DE-AC02-76CH03000) and 



by NASA grant NAG5-10842. D.H was supported by 
DOE grant DE-FG02-95ER40896, J.G.L. and S.P. by 
DOE grant DE-FG03-94ER40833 and T.J.W. by DOE 
grant DE-FG05-85ER40226. 



[1] 



[2] 
[3] 

[4] 

[5] 
[6] 
[7] 
[8] 
[9] 



[10] 

[11] 

[12] 
[13] 

[14] 

[15] 

[16] 
[17] 



[18] 
[19] 
[20] 

[21] 



L. Wolfenstein, Nucl. Phys. B 186, 147 (1981); S. T. Pet- 
cov, Phys. Lett. B 110, 245 (1982); S. M. Bilenky 
and B. Pontecorvo, Sov. J. Nucl. Phys. 38, 248 (1983); 
S. M. Bilenky and S. T. Petcov, Rev. Mod. Phys. 59, 
671 (1987); C. Giunti, C. W. Kim and U. W. Lee, Phys. 
Rev. D 46, 3034 (1992); J. P. Bowes and R. R. Volkas, 
J. Phys. G 24, 1249 (1998); K. R. Balaji, A. Kalliomaki 
and J. Maalampi, Phys. Lett. B 524, 153 (2002). 
M. Kobayashi and C. S. Lim, Phys. Rev. D 64, 013003 

(2001) . 

R. Foot, H. Lew and R. R. Volkas, Phys. Lett. B 272, 67 
(1991); Mod. Phys. Lett. A 7, 2567 (1992); R. Foot and 
R. R. Volkas, Phys. Rev. D 52, 6595 (1995). 
V. Berezinsky, M. Narayan and F. Vissani, Nucl. Phys. 
B 658, 254 (2003). 

A. Aguilar et al, Phys. Rev. D 64, 112007 (2001). 
P. Langacker, Phys. Rev. D 58, 093017 (1998). 
G. Cleaver et al., Phys. Rev. D 57, 2701 (1998). 
D. Chang and O. C. Kong, Phys. Lett. B 477, 416 (2000). 
K. R. Dienes, E. Dudas and T. Gherghetta, Nucl. Phys. 
B 557, 25 (1999); N. Arkani-Hamed, S. Dimopoulos, 
G. R. Dvah and J. March-Russell, Phys. Rev. D 65, 
024032 (2002); G. R. Dvali and A. Y. Smirnov, Nucl. 
Phys. B 563, 63 (1999). 

T. J. Weiler, W. A. Simmons, S. Pakvasa and 
J. G. Learned, hep -ph/94114 32 D. J. Wagner and 
T. J. Weiler, Mod. Phys. Lett. A 12 (1997) 2497. 
R. M. Crocker, F. Melia and R. R. Volkas, Astrophys. J. 
Suppl. 130, 339 (2000); Astrophys. J. Suppl. 141, 147 

(2002) . 

C. Lunardini and A. Y. Smirnov, Nucl. Phys. B 583, 260 
(2000); Phys. Rev. D 64, 073006 (2001). 
A. D. Dolgov et al, Nucl. Phys. B 632, 363 (2002); 
K. N. Abazajian, J. F. Beacom and N. F. Bell, Phys. 
Rev. D 66, 013008 (2002); Y. Y. Wong, Phys. Rev. D 
66, 025015 (2002). 

J. G. Learned and S. Pakvasa, Astropart. Phys. 3, 267 

(1995); L. Bento, P. Keranen and J. Maalampi, Phys. 

Lett. B 476, 205 (2000); H. Athar, M. Jezabek and O. Ya- 

suda, Phys. Rev. D 62, 103007 (2000). 

J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa 

and T. J. Weiler, Phys. Rev. Lett. 90, 181301 (2003); 

hep-ph/0309267 [Phys. Rev. D (to be published)]. 

J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa and 

T. J. Weiler, Phys. Rev. D 68, 093005 (2003). 

S. N. Ahmed et al., nucl-ex/0309004 K. Nishikawa, 

talk at Lepton-Photon 2003 Conference, August 2003, 

http: / /conferences. fnal.gov /lp2003 / 

W. L. Freedman et al., Astrophys. J. 553, 47 (2001). 

L. Stodolsky, Phys. Lett. B 473, 61 (2000). 

S. Choubey and S. F. King, Phys. Rev. D 67, 073005 

(2003) . 

P. Keranen, J. Maalampi, M. Myyrylainen and J. Riitti- 
nen, hep-ph/0307041j 



